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We present a quantitative many-body analysis using the GW approximation of the decay rate F 
due to electron-electron scattering of excitations in the Shockley surface state band of Ag(lll), as 
measured using the scanning tunnelling microscope (STM). The calculations include the perturbing 
influence of the STM, which causes a Stark-shift of the surface state energy E and concomitant 
increase in F. We find F varies more rapidly with E than recently found for image potential states, 
where the STM has been shown to significantly affect measured lifetimes. For the Shockley states, 
the Stark-shifts that occur under normal tunnelling conditions are relatively small and previous 
STM-derived lifetimes need not be corrected. 
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The femtosecond lifetimes r of electronic excitations 
in noble metal surface bands"'^ can be determined us- 
ing the scanning tunnelling microscope (STM). The 
STM-based techniques complement those that use photo- 
excitation such as photoemission spectroscopy (PES) and 
two-photon photoemission^, being applicable to both 
electron-like and hole-like excitations in combination 
with the advantages of atomic-scale spatial resolution. 
There are two principal methods available. The first is 
based upon line shape analysis of differential conductivity 
(dl/dV) measurements^!^, along with artificial or natu- 
rally occurring nanoscale resonator structures which are 
used to induce spectral structure over a wide range of 
energies^i^. The second approach is based upon the 
spatial variation of quantum interference patterns that 
are visible in dl / dV measurementsS'S*^. These are anal- 
ysed to determine the phase coherence length L^, which 
can be converted to the lifetime via the group velocity 

Vg-. T = L^/Vg. 

An analysis of Ag(lll) Shockley lifetime values de- 
termined using these two different STM techniques has 
demonstrated that considerable agreement exists, over 
a range of energies^^. The lifetimes are also compara- 
ble to state-of-the-art many-body calculations using the 
GW method, but decrease more rapidly with increasing 
energy. In comparing theoretical and experimental life- 
times no consideration has yet been given to the possible 
consequences of the perturbing influence of the STM tip, 
which is known to cause a measurable Stark-shift in the 
surface state energiesii. Many-body GW calculations 
have shown that in the case of the higher-lying image 
potential surface states this Stark-shift is accompanied 
by a significant increase in the inelastic decay ratei^i; the 
electric field between the STM tip and sample leads to a 
doubling of the decay rate under normal tunnelling con- 
ditions. Here we perform a quantitative study using the 
many-body GW method to quantify for the first time 
the impact of the STM tip on the electron-electron scat- 
tering rate of Shockley state electrons. We discount the 



possibility of a significant Stark-shift induced change to 
the electron-phonon contribution, which has previously 
been shown to be constant for excitations with energies 
in excess of 20 meVii^ 

Our calculations are based upon the approach de- 
veloped by Chulkov and coworkergi^, and used widely 
in calculations of surface state dynamics^ including the 
lifetimes of Stark-shifted image potential statesii^ The 
damping rate or inverse lifetime of an excitation in the 
state ip{r) with energy E is obtained from the expec- 
tation value of the imaginary part of the electron self- 
energy, T,{r,r';E): 



r = r-i = 



= -2 JdrJdr'i;*{r)lmJ:{r,r';E)i(j{r'). (1) 

In the GW approximation— of many-body theory the 
imaginary part of the self energy is calculated in terms 
of the screened Coulomb interaction W and the Green 
function G 

ImE(r,r';e) = / rfe' ImG(r, r'; e')Imiy(r, r'; e-e'). 

71" Jef 

(2) 

A successful account of the decay rates of the noble metal 
Shockley surface states is possible using for the Green 
function its non-interacting counterpart 



G{r,r';e)^J2 



e-E,+ iQ+ ' 



(3) 



and evaluating the screened Coulomb interaction in the 
random phase approximation (RPA), 



W{r,r'\Lo) = v{r — r') + J dri J dr2v(r — ri) 

xx°{ri,r2;cj)Wir2,r';Lj) (4) 

where v is the bare Coulomb interaction and x° is the 
density-density response function of the non-interacting 
electron system. 
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FIG. 1: Surface band structure of Ag(lll). The shaded 
area and solid line are the projected continuum and surface 
state dispersion from ab-initio calculations which we use for 
calculating the lifetime. Symbols are STM-derived dispersion 
data: circles, Ref. l2fll : squares, Ref. j HI The dashed lines are 
parabolic dispersions with effective masses m* = 0.42 (surface 
state band) and 0.24 (band edge). 



For the evaluation of the single particle states we have 
solved the Schrodinger equation using standard pseu- 
dopotentials varying in the direction perpendicular to the 
surfaceiSiii. The pseudopotential does not describe the 
d-electrons of the substrate, which are anyway too low 
in energy to play a significant role as final states for the 
decay, but their contribution to the screening is included 
via a polarisable backgroundiSiiS. 

Previous calculations of the decay rate of the Ag(lll) 
Shockley state have assumed parabolic dispersion with 
effective masses of m* = 0.42 and m* = 0.24 for the 
intrinsic surface state and the lower band edge respec- 
tively, but we find that over the extended energy range 
of interest here these provide a poor description of the 
dispersion that is found in ab-initio calculations, as may 
be seen in Fig. ^ The parabolic dispersions result in the 
surface state crossing the band edge, which is incorrect. 
In the calculations reported here we have therefore used 
the ab-initio dispersions. As previously notedSLSS there 
are important changes in the shape of the surface state 
wave function with k, the surface wave vector, and we 
take these into account by recalculating the wave func- 
tion for different k with the pseudopotential parameters 
changed to take into account the appropriate ab-initio 
band edges and surface state energy. Finally we have 
also modified the potential outside the surface in order 
to take account both of the multiple images present in 
the tunnel junction geometry and the electric field be- 
tween tip and sample due to the applied bias voltage, as 



described in Ref. ■ There it was found that measured 
Stark-shifts were best described by assuming no contact 
potential existed between tip and sample, and so in the 
present calculations we have assumed that the tip and 
sample have similar work function. 

Wc first consider the influence of the STM tip on the 
lifetimes of holes at the bottom of the Ag(lll) surface 
state band (fc — 0), which have been determined exper- 
imentally from line shape analysis of the step-like onset 
that is seen in open- feedback dl/dV measurements on de- 
fect free regions of the surfaced The onset occurs when 
the bias voltage coincides with the surface state binding 
energy, so that the Fermi energy of the tip coincides with 
the minimum of the surface state band. This is Stark- 
shifted downwards from the value of Eq — —63 meV ob- 
served in PES experiments^ by the STM-derived electric 
field. The field itself varies according to the tip-sample 
separation, which depends experimentally upon the size 
of the current prior to opening the feedback loop, so that 
different choices for this current result in different shifts. 
Limot et ali^ were able to follow the change in Eo from 
—66 to —80 meV using currents in the range 50 pA — 6 
^A. 

We simulate these conditions as follows. For a given 
tip-sample separation we solve the Schrodinger equation 
for the tunnel junction to determine the energy of the 
surface state, varying the bias voltage until the calcu- 
lated Stark-shifted binding energy and the bias voltage 
coincide (threshold tunnelling condition). At this point 
we use equations ©-(0J to calculate the hole decay rate 
r that applies for this value of the Stark shift, before re- 
peating for different values of the tip-sample separation. 

The results are shown in Fig. [21 Our field-free value 
agrees with the 2.8 meV lifetime broadening calculated 
by Garci'a-Lekue et ali^ using similar methods. Includ- 
ing the field due to the STM tip, we find a linear increase 
in the decay rate as the surface state is Stark-shifted to 
greater binding energies, with F some 30% greater when 
the binding energy is Eq = —80 meV than at the field- 
free energy of —63 nieV. The rate of increase is com- 
parable to that found for image state electrons, with 
\dr/dEo\ = 0.050, compared to dT/dEi = 0.037 for the 
Cu(OOl) n = 1 image statCf-'^^ although the range of ac- 
cessible Stark-shifted energies is significantly smaller so 
that the absolute change in decay rates that occurs for 
the Shockley state is much smaller. The integral in Eqn. 
^ is over final states, and by isolating contributions to 
the imaginary self energy due to the surface state band 
and using just those to calculate the decay rate allows F 
to be decomposed into intraband and interband contribu- 
tions. We find that the variation in the decay rate shown 
in Fig. 121 is primarily due to the change in the rate of in- 
traband scattering, the process by which the hole is filled 
by an electron from within the surface state band itself. 
The interband contribution, in which the hole is filled 
by an electron from the bulk continuum, contributes less 
than 6% of the overall decay rate. 

To understand the origin of the change in the decay 
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FIG. 2: Decay rate F due to electron-electron scattering of 
hole-like excitations at the bottom of the Shockley surface 
state band on Ag(lll), as a function of Stark-shifted energy 
Eq. Filled triangles show the total decay rate, and open tri- 
angles the intraband contribution. The inset illustrates intra- 
band and interband decay processes. 
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FIG. 3: Calculated change in the decay rate induced by the 
electric field of the STM, for Ag(lll) surface state electrons 
at the tunnelling threshold {E — Ef = eV). A tip-sample 
separation of 8 A has been used. The inset shows the size of 
the Stark-shift. 



rate we have performed a number of auxiliary calcula- 
tions, focusing on the dominant intraband decay rate. 
These show that the physical basis for the increase in the 
decay rate is the increased number of final states available 
for decay. The calculations are as follows. At the field- 
free surface we find Tintra — 2.66 meV, and when the 
surface state is Stark-shifted to a binding energy of —80 
meV we find Tintra = 3.41 meV. If we calculate the decay 
rate using the Stark-shifted eigenvalues but replace the 
corresponding wave functions by their unperturbed coun- 
terparts, we find Tintra = 3.37 meV, which is close to the 
result of the calculation using the field-perturbed wave 
functions, and conversely using the Shockley wave func- 
tions calculated in the presence of the applied field, but 
the unperturbed energies, gives Tintra = 2.67 meV, close 
to the field-free decay rate. Hence changes in the wave 
functions of the surface and bulk states due to the STM- 
induced electric field are of minor importance, which is 
different to the situation found for the n — 1 image state 
at Cu(OOl) where changes in the wave function are the 
most important contributor to the increased decay raters- 
Next we have also calculated the decay rate using the 
field-free energies and wave functions, but with the sur- 
face state wave function in the Green function Eqn. Q 
scaled by a/80/63; 80/63 is the ratio of the number 
of electrons in the surface state band in the perturbed 
and unperturbed systems. The resulting decay rate is 
Tintra = 3.41 mcV. This simple re-normalisation of the 
contribution of the surface state to the sum (integral) 
over final states in the evaluation of the self-energy, Eqn. 



0, is sufficient to reproduce the result of the calcula- 
tion using the properly Stark-shifted states, and indi- 
cates that it is only the change in the total number of 
states available to fill the hole that matters, and not their 
changed energies. 

We have extended this study to also consider the de- 
cay rate of excitations at energies further up the Ag(lll) 
Shockley surface state band. These have been measured 
experimentally using line shape analysis in nanoscale 
resonators^ and from studies of quantum interference 
patterns near steps^ and in triangular adatom corralsi^ 
We consider energies between and 1 eV where the sur- 
face state is well defined. Physically this energy range is 
of interest as within it there is a cross-over in the primary 
decay channel from intraband to interband dominated 
decayj^ In these calculations we keep the tip-sample sep- 
aration fixed at 8 A as the bias voltage is varied (this 
distance results in a 4 meV Stark-shift of the band edge 
state, as was presented in Ref. 3]), and we once again 
calculate the decay rate of the surface state at the tun- 
nelling threshold, where the energy of the Stark-shifted 
state coincides with the Fermi energy of the tip. 

The results are presented in Fig. |3| Using this fixed 
tip-sample separation we find that the Stark-shift AE of 
the surface state energy increases linearly with the bias 
over the range of voltages considered, and reaches 22 meV 
when the bias is 1 V (Fig.O inset However, accompa- 
nying the field-induced modification of the surface state 
we find a rather small change in the decay rate. The 
additional numerical complexity associated with the re- 
duced symmetry of an excitation with a finite wave vector 
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results in a greater uncertainty in the calculated decay 
rate away from the band edge. Allowing for this, we find 
that the change in the decay rate is less than 0.3 nieV 
over most of the range, only reaching 0.7 ± 0.15 mcV for 
a bias of 1 V. 

We now consider our results as a whole and in the con- 
text of experimental studies. In the case of holes at the 
bottom of the surface state band, in the measurements 
of Kliewer et al^ the onset was observed at —67 meV for 
Ag(lll), corresponding to a Stark-shift of 4 meV from 
the PES field-free valueiS^ Our calculations reported here 
(Fig. [2Jl show that the the corresponding decay rate is 
increased by 0.2 meV above the field-free value, which is 
within the experimental uncertainty^ For electron-like 
excitations in the surface state band the decay rate in- 
creases with energy due to the increase in available decay 
channels, and for an energy 1 eV above Ep the corre- 
sponding contribution to the linewidth is approximately 
40 meV. The STM- induced increase of 2± 0.7 nieV is small 
in comparison, and again well within the experimental 
uncertaintyiiSi 

Therefore, to conclude, on the basis of our quantitative 
theoretical many-body analysis using the GW approx- 



imation we have demonstrated that in contrast to the 
case of image potential states, the lifetimes of excitations 
in noble metal Shockley surface state bands determined 
using the STM are not sufficiently affected by the electric 
field of the STM tip that previously reported values need 
to be corrected. We have identified that relative changes 
of up to 30% can occur under extreme but experimen- 
tally accessible tunnelling conditions, and the decay rate 
actually changes more rapidly as a function of the Stark- 
shift for Shockley states than for image potential states. 
However the significantly smaller Stark-shifts that occur 
when tunnelling via the lower-lying Shockley states un- 
der normal tunnelling conditions mean that the absolute 
changes in the decay rate lie within experimental uncer- 
tainties previously reported. 
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